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ABSTRACT 



Context. Usually, important parameters of young, low-mass star clusters are very difficult to obtain by means of pho- 
tometry, especially when differential reddening and/or binaries occur in large amounts. 

Aims. We present a semi-analytical approach (ASA m i n ) that, applied to the Hess diagram of a young star cluster, is 
able to retrieve the values of mass, age, star-formation spread, distance modulus, foreground and differential reddening, 
and binary fraction. 

Methods. The global optimisation method known as adaptive simulated annealing (ASA) is used to minimise the residu- 
als between the observed and simulated Hess diagrams of a star cluster. The simulations are realistic and take the most 
relevant parameters of young clusters into account. Important features of the simulations are: a normal (Gaussian) dif- 
ferential reddening distribution, a time-decreasing star-formation rate, the unresolved binaries, and the smearing effect 
produced by photometric uncertainties on Hess diagrams. Free parameters are: cluster mass, age, distance modulus, 
star-formation spread, foreground and differential reddening, and binary fraction. 

Results. Tests with model clusters built with parameters spanning a broad range of values show that ASA m i„ retrieves 
the input values with a high precision for cluster mass, distance modulus and foreground reddening, but somewhat 
lower for the remaining parameters. Given the statistical nature of the simulations, several runs should be performed 
to obtain significant convergence patterns. Specifically, we find that the retrieved (absolute minimum) parameters con- 
verge to mean values with a low dispersion as the Hess residuals decrease. When applied to actual young clusters, the 
retrieved parameters follow convergence patterns similar to the models. We show how the stochasticity associated with 
the early phases may affect the results, especially in low-mass clusters. This effect can be minimised by averaging out 
several twin clusters in the simulated Hess diagrams. 

Conclusions. Even for low-mass star clusters, ASAmin is sensitive to the values of cluster mass, age, distance modulus, 
star-formation spread, foreground and differential reddening and, to a lesser degree, binary fraction. Compared with 
simpler approaches, the inclusion of binaries, a decaying star-formation rate and a normally distributed differential 
reddening, appear to yield more constrained parameters, especially the mass, age and distance from the Sun. A robust 
determination of cluster parameters may have a positive impact on many fields. For instance, age, mass and binary 
fraction are important for establishing the dynamical state of a cluster, or deriving a more precise star-formation rate 
in the Galaxy. 

Key words. (Galaxy:) open clusters and associations; Galaxy: structure 
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Several events occurring on short timescales combine to 
make the first few 10 7 yr the most critical period in a star 
cluster's life, especially for the poorly-populated embedded 
clusters (ECs). Driven mainly by the impulsive removal of 
the parental molecular gas by supernovae and massive-star 
winds, the rapid and deep changes in the cluster internal 
dynamics lead to the escape of varying fractions of member 
stars to the field. The end result is that most of the low- 
mass clusters are dis solved before reaching ~ 40 Myr of age 



into gravitationally bound o pen clusters (e.g. |Lada fc Lada 
2003| |Bonatto fc Bica||2011| ) 



Low-mass clusters undergoing such a rapidly chang- 
ing phase are usually characterised by an under-populated 
main sequence (MS) and significant fractions of pre-main 
sequence (PMS) stars, all mixed up with a spatially non- 
uniform dust and gas distribution. To complicate mat- 
ters, star formation within a cluster appears to extend 
for a period that may be comparable to the cluster age 



(e.g. Stauffer et al. 19971, while the MS and PMS stars 



(e.g. Tutukov 1978 



Goodwin fc Bastian||2006| . This sce- 
nario is consistent with current estimates suggesting that 
less than ~ 5% of the Galactic ECs dynamically evolve 
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may be a rranged in (unk nown) fractio ns of binaries (e.g 



Weidner, Kroupa & Maschberger 20091. Besides, the vast 
majority of the binary pairs cannot be resolved in Colour- 
Magnitude Diagrams (CMDs). As a result, the CMDs 
of young, low-mass star clusters tend to be very diffi- 
cult to interpret. And, consequently, the determination of 
fundamental cluster parameters (e.g. age, distance, mass, 
foreground and internal reddening, star-formation spread 



1 



C. Bonatto et al.: Hidden properties of young star clusters 



(SFS), binary fraction, etc) may be somewhat subjec- 
tive and unreliable, particularly when only photometry 
is available to work with. Examples of clusters evolving 
along the early phase and characterised by such complex 
CMDs are abundant in the recent literature. F or instance: 
NGC6611 dBonatto Santos Jr. k Bica||2006|) , NGC4755 



( Bonatto et al.||2006D, NGC2244 QBonatto fc Bic a 2009a), 
Bochum l (|Bica, Bonatto k Dutra| |2008"| , Pismis5 and 
vdB80 dBonatto k Bica| |2009b[), CollinSer 197 and vdB 92 
( Bonatto k Bica||2010a| ) 

More broadly, a robust determination of fundamental 
parameters of young clusters may provide important con- 
straints for studies dealing with dynamical state and clus 



ter dissolution time-scales (e.g. Goodwin 2009 



Baumgardt k Giclcs 2010) , infant m ortality (e.g 



Lamcrs 
Lad a k 



Lada 2003; Goodwin k Bastian 2006 1 , star- formation rate 



(SFR) in the Galaxy (e.g. Lamers k Gieles 2006 Bonatto 



k Bica|[20T7T ), among others. 



Over the years, several approaches have been proposed 
to tackle the important task of fi nding reliable parameters 
of young clusters. Among these, |Naylor k Jeffries| ([2006) 
employ a maximum-likelihood method to Hess diagram^ 
simulations (including binaries) to derive distances and 
ages. However, they do not consider differential reddening 
(DR), and their method appears to be more efficient for 
clusters older than ~ 30 Myr. Hillenbrand, Bauermeister 



k White ( 2008 ) model CMDs of star forming regions and 
young star clusters by means of varying star-formation his- 
tories; because of confusion between signal and noise in 
CMDs, onl y marginal evidence for moderate ag e spreads 
was found. |da Rio, Gouliermis k Gennaro (|2010 ) ad d PR , 
age spreads and PMS stars to the Naylor k Jeffries (2006) 
method, but adopt distance and redd ening values from pre 
work 



vious 

to estimate 



|Stea d 
the" 



Hoare (2011) employ Monte Carlo 



age of ECs observed with near-infrared 
(UKIDSS) photometry. More recently, Bonatto, Bica k 
Lima (2011) (hereafter Paper I) included the smearing ef- 
fect of photometric uncertainties on CMDs to approach the 
above issue by means of brute-force. However, because of 
computer limitations, the binary fraction, SFS, SFR, and 
the shape of the DR distribution had to be taken as fixed 
parameters. 

In summary, the problem of obtaining reliable cluster 
parameters clearly lacks a comprehensive method that in- 
cludes most (if not all) of the relevant conditions and pa- 
rameters usually associated with the early cluster phases. 
Here we present a semi-analytical approach that follows 
along this direction. Instead of relying on brute-force, we 
now make use of the global optimisation method kno wn as 
adaptive simulated an nealing (ASA 



adapted from Goffe 
Ferrier k Rogers] |T994) to search for the set of parameters 



that best reproduces the photometric properties (i.e., the 
CMD or Hess diagram) of a cluster. The free parameters 
are the apparent distance modulus, foreground reddening, 
cluster mass, age, SFS, binary fraction, and DR, which are 
used to build the simulated Hess diagram for a given cluster. 
The function to be minimised by ASA is related to differ- 
ences between the observed and simulated Hess diagrams. 
For conciseness, hereafter we will refer to our approach as 



1 Hess diagrams contain the relative density of occurrence 
of stars in different colour-magnitude cells of the Hertzsprung- 




Russell diagram ( pless|[l924| 



t SF (Myr) 

Fig. 1. Top: The SFR (mass converted into stars with re- 
spect to the cluster mass) as a function of star-formation 
time for a cluster undergoing continuous star formation for 
20 Myr. Bottom: Same as above for the number of formed 
stars. 

By explicitly including the binary fraction, SFS and the 
shape of the DR distribution as free parameters on a semi- 
analytical approach, this work supersedes the procedure 
highlighted in Paper I. 

The present paper is organised as follows. In Sect. [2] we 
discuss the relevant effects that introduce observational lim- 
itations on CMDs. In Sect. [3] we provide details on the clus- 
ter simulations. In Sect. [4] we briefly describe the ASA min 
minimisation method, test it on model clusters and discuss 
the strategy adopted to obtain robust parameters. In Sect. [5] 
we apply ASA min on actual young clusters, contextualise 
the results and discuss the low-mass cluster stochasticity. 
Concluding remarks are given in Sect. [6] 

2. Intrinsic limitations of CMDs 

Most of the problems related to obtaining reliable funda- 
mental parameters of young clusters are illustrated in Fig.[T] 
of Paper I, in which a model CMD is subject to different val- 
ues of distance modulus, DR and binary fraction. Clearly, 
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when a moderate value of DR is added to the model CMD, 
the scatter and spread in the stellar sequences (resulting 
also from photometric uncertainties and binaries) tends to 
mask any intrinsic relationship between parameters of in- 
dividual stars (e.g. the age t*, mass m*, magnitudes and 
colours) and those of the isochrones, an effect that increases 
significantly with distance modulus, DR amount and binary 
fraction. 

Following Paper I, we keep working with 2MAS^] pho- 
tometry and, in particular, with J x ( J — K s ) CMDs (the 
colour least affected by photo metric errors and the best dis- 
criminant for PMS stars - e.g. Bonatto fc Bica|[2010b[ ) and 
the corresponding Hess diagrams. Reddening transforma- 
tions are based on the absorption relat ions Aj/Ay = 0.276 
and AkJAv = 0.118, with R v = 3.1 flCardelli, Clayton fc 



|Mathis||1989| |Dutra, Santiago fc Bica||2002p . To avoid con- 
fusion between notations, we refer to the foreground red- 
dening as E(J ~ K s ), where E(J — K s ) = 0.158A V , and to 
the differential reddening as SAy. 

As shown in Paper I, the photometric uncertainties play 
an important role in shaping the CMD (or Hess diagram) 
morphology. Here we take them explicitly into account 
by assuming a normal (Gaussian) distribution of errors. 
Formally, if the magnitude (or colour) of a star is given by 
x ± a x , the probability of finding it at a specific value x is 

given by P(x) — — e~ 5 (^^) . Thus, for each star we 

compute the fraction of the magnitude and colour that falls 
within a given cell of the Hess diagram (i.e., the cell den- 
sity). Essentially, this corresponds to the difference of the 
error functions computed at the cell borders. By definition, 
the sum of the colour and magnitude density over all Hess 
cells should be the number of input stars. As a compromise 
between CMD resolution and computational time, the Hess 
diagrams used in this work consist of magnitude and colour 
cells of size AJ = 0.2 and A( J — K s ) = 0.02, respectively. 

As a photometric quality criterion, we only keep stars 
with errors lower than 0.2, which also serves to emulate 
the photometric completeness function of the observations. 
The effect of this restriction on the luminosity function of 
a model cluster containing an arbitrarily large (for statis- 
tical purposes) number of stars can be seen in Fig. [T] of 
Paper I. After following the steps described below for as- 
signing the mass, magnitudes, colours and uncertainties to 
each star, we apply the error restriction of keeping only 
stars with errors < 0.1 and 0.2. Compared to the complete 
luminosity function, which presents a turnover at J 16.9, 
the restricted functions have turnovers that smoothly shift 
towards brighter magnitudes. 



3. A 



simulation recipe for young clusters 



Major changes with respect to the procedures adopted in 
Paper I for simulating young clusters are as follows: (i) The 
DR is described as a normal distribution, characterised by 
the mean SAy and standard deviation around the mean 
cr(SAy), where both are free parameters, (ii) The SFR is a 
linearly decreasing function of time, which is a more realis- 
tic a ssumption, especially for cluste rs older than ~ 10 Myr 
(e.g. Belloche et al. 2011 Bate 2011 and references therein). 
Although somewhat arbitrary, this assumption has the ad- 
vantage of requiring a single free parameter. In this case, it 



2 The Two Micron All Sky Survey, All Sky data release 
( Skrutskie et al.|2006 l 



is just the cluster age (t c i u ), which we assume as the begin- 
ning of the star-forming event. For a cluster of mass M c i u 
and age t c i u , the SFR is the mass M converted in stars as 
a function of time after the star formation start t s f , which 
we express as |^ = (1 - t sf lt clu ). This SFR is il- 

a isf i^clu 

lustrated in Fig. n] (top panel) , for a 20 Myr old cluster that 
has been forming stars during the same time. Equivalently, 
the fraction of formed stars is also a linearly decreasing 
function of time (bottom), (in) The SFS timescale (r s / s ) is 
now a free parameter; instead of a fixed (and assumed to be 
equal to the cluster age) value, it can vary within the range 
< T s f s < t c i u . Specifically, t s / s corresponds to the differ- 
ence between t c i u and the minimum (isochrone) age com- 
patible with the CMD stellar distribution, (iv) Binaries are 
expected to survive the early evolutionary phase of low- 
mass clusters, with the unresolved pairs being somewhat 
brighter and redder than the single stars, thus resulting 
in some broadening of the CMD stellar seq uences, espe- 
cially the MS (e.g. |Naylor fc Jeffries] |2006] ). We include 
them in our simulations by means of the (free) parame- 
ter fbin, which measures the fraction of unresolved bina- 
ries in a CMD. According to this definition, a CMD with 
Nqmd detections and characterised by the binary fraction 
fbin, would have a number of individual stars expressed as 
iV* = (1 + fbin)NcMD- As an additional assumption, bi- 
naries are formed by pairing stars with the closest ages, 
regardless of the individual masses. This gives rise to a sec- 
ondary to primary stellar mass-ratio (q = m s /m p ) that 
smoothly increases from very-low values up to q ~ 0.25, 
and decreases for higher values of q (Fig. [I] of Paper I) . 

Thus, to simulate the CMD of a cluster on a given 
photometric system we: (i) Start with an artificial clus- 
ter of mass M c i u , age t c i u , SFS timescale r s / s , appar- 
ent distance modulus (m — M)j, foreground reddening 
E( J — K s ), mean DR SAy and standard deviation a(SAy), 
and binary fraction fbin- (H) Assign each star an age accord- 
ing to t c i u — r s f s < t* < tdu, with a probability following a 
linearly decreasing age distribution (Fig. [I]). (Hi) Select the 
indiv idual ste l lar m asses (to*) by randomly taking values 
from Kroupa (2001) mass function (with the upper mass 
value consistent with the age), until J2i m *i = ^clu- (iv) 
Compute the magnitudes and colours of each star according 
to t* and to* , and the mass to light relations taken from the 
isochrone corresponding to i*. (v) Form the (unresolved) 
binary pairs according to fu n and compute their magni- 
tudes and colours, (vi) For each star, randomly compute 
DR from a normal distribution characterised by the mean 
SAy and standard deviation a (SAy). (vii) Apply shifts in 
colour and magnitude according to the values of E( J — K s ) 
and (m — M)j. (viii) Assign each artificial star a photo- 
metric uncertainty based on the average 2MASS errors and 
magnitude relationship, (ix) For more realistic representa- 
tiveness, add some photometric noise to the stars. This step 
is taken to minimise the probability of stars with the same 
mass having exactly the same observed magnitude, colour, 
and uncertainty in the CMD. To do this, consider a star 
(of mass to* and age t*) with an intrinsic (i.e., measured 
from the corresponding isochrone) magnitude mag and as- 
signed uncertainty a mag . The noise-added magnitude mag 
is then randomly computed from a normal distribution with 
a mean mag and standard deviation o~ mag . (x) Apply the 
same detection limit to the model CMD as for the observa- 
tions, so that model and data share a similar photometric 
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completeness function; in practise, this means that stars 
with photometric errors higher than 0.2 are discarded, (xi) 
Finally, build the corresponding Jx ( J — K s ) Hess diagram. 

To minimise the critical sto chast icity associated with 
low-mass clusters (see, e.g. Sect. 5.1), steps (ii) to (x) are 
repeated for N S i m (twin) clusters of mass M c i u , age t c i Ul 
SFS timescale r s / s , and binary fraction fun- Thus, each 
cell of the artificial Hess diagram (step (xi)) contains the 
density of stars averaged over the N sim clusters. Some tech- 
nical details are described below. 



The broken mass function of Kroupa (2001) is de- 
fined as dN/dm oc mT^ 1+x \ with the slopes \ — 0-3 
for 0.08 < m(M Q ) < 0.5 and \ = 1-3 for m(M ) > 
0.5. The relations of photometric uncertainty with mag- 
nitude for the J and K s 2MASS bands are well repre- 
sented by a j = 0.0214 + 2.48 x 10" 8 exp (J/1.071) and 
a Ks = 0.0193 + 9.59 x 10~ 9 exp (1^/1. 067). These relations 
apply to the range 3 < J, K s < 20. 

The stellar mass/luminosity relation is t aken from the 
solar- met allicity isochrone sets of Padovcp ] (Girardi et al 



2002) and |Siess Dufour k Forestini| pJOOf T Both sets have 
been merged, since Padova isochrones should be used only 
for the MS (or more evolved sequences) , and those of Siess 
apply to the PMS^j The merging occurs at the MS entry 
point, at 6.5 M for the isochrones younger than 8Myr, 
5.5 M for lOMyr, 4.5 M Q for 20Myr, and 3.5 M for 
> 30 Myr. We now work with a rather high time resolu- 
tion, considering isochrones with ages from to 10 Myr 
(with a step of lMyr), 10 to 20 Myr (step of 2 Myr), and 
20 to 50 Myr (step of 5 Myr) . The cluster models consist of 
stars more massive than 0.1 Mq, the lowest available mass 
in the PMS isochrones. Magnitudes, colours and mass for 
model stars with intermediate age values are obtained by 
interpolation among the neighbouring isochrones. Also, the 
maximum stellar mass present in the adopted isochrones 
ranges from 60 M (at 0.2 Myr) through 36 M (5 Myr), 
19 M (10 Myr), 9M (30 Myr), and 7.3 M (50 Myr). 



4. Parameter optimisation with ASA m i n 

As discussed in the previous section, the problem 
now is to search for the set of 8 free parameters P = 
{M c i u ,t clu , T sfs ,(m- M)j,E{J ~ K s ),lA v ,cr(6A v )J bm } 
that produces the best match between the observed and 
simulated Hess diagrams of a given star cluster. This is 
achieved by minimising the root mean squared residual 
between both diagrams (i? rms ) that, by construction, is 
a function of the free parameters, i.e., R rm s — Rrms{P)- 
As shown in Paper I, such a task would take an extremely 
long time for the brute-force method, especially when one 
considers a realistic (relatively wide) range of values for 
such an extended set of parameters. 

Among the several optimisation methods available in 
the literature, the adaptive simulated annealing (ASA) ap- 
pears to suit our purp oses, because it is relatively ti me ef- 
ficient and robust (e.g. Goffe, Ferrier fc Rogers|1994 ). ASA 



3 Built for the 2MASS filters at http://stev.oapd.inaf.it/cgi- 
\bin/cmd^ 

* As a caveat, some colour bias may be present in the Siess 
isochrones, because they u se a single T e -colour relation (from 



has been shown to be a global optimisation technique that 
distinguishes between different local minima, and the R rms 
hyper-surface is characterised by the presence of such fea- 
tures (Paper I). Simulated annealing derives its name from 
the metallurgical process by which the controlled heating 
and cooling of a material is used to increase the size of its 
crystals and reduce their defects. If an atom is stuck to 
a local minimum of the internal energy, heating forces it 
to randomly wander through higher-energy states. In the 
present context, a state is a simulation corresponding to 
a specific set of values of the parameters to be optimised. 
Then, the slow cooling increases the probability of finding 
states of lower energy than the initial one. 

Minimisation starts by the definition of individual 
search-ranges (wide enough to permit the occurrence of 
any value compatible with the cluster nature) and varia- 
tion steps for each free parameter; next, an initial point is 
randomly selected and the starting residual value R l rms is 
computed. At this point, we define the initial temperature 
as T = R l rms - Then ASA takes a step (i.e, changing the 
initial parameters) and the new value -R^ms i s evaluated. 
Specifically, this implies that a new Hess diagram has to 
be simulated with the new parameters. By definition, any 
downhill (-ffj-ros < Rrms) ste P i s accepted, with the process 
repeating from this new point. However, uphill moves may 
also be taken, with the decision made by the Metropolis 



(Metropolis et al. 1953) criterion, which allows small up- 



hili moves while rejecting large ones, thus enabling ASA 
to escape from local minima. When an uphill move is re- 
quired, ASA computes the value of V = e~( A /' 1 ~\ where 
A = Jij-ms — Rrms- A/T is positive in an uphill move and 
so, V is a number between and 1 that is compared with 
a random number < M < 1. If V > A/", the uphill move 
is accepted and the algorithm moves on from that point; 
in case of rejection, another point is chosen for a trial eval- 
uation of V and M . Clearly, large values of A and low T 
make acceptance of an uphill move less likely. After each 
successful move, the temperature of the system is reduced 
according to T — > f 7", which makes uphill moves less likely 
to be accepted as ASA focuses upon the most promising 
area for optimisation. For a more precise (but slower) con- 
vergence rate, we use / = 0.95. Variation steps decrease 
as the minimisation is successful and ASA closes in on the 
global minimum (P —¥ P m i n ). The termination criterion 
for a run occurs when A < 10~ 6 or, for runtime sake, the 
number of trial function evaluations - for a single uphill 
move - reaches 2.5 x 10 5 . 

We define the root mean squared residual between ob- 
served (H b s ) and simulated (H S i m ) Hess diagrams, com- 
posed of n c and n m colour and magnitude cells, as 



R, 



\ 



1 [H obs (i,j) 



E 



H sim {i,j)Y 



Nobs .4< [H obs (i, j) + H sim (i,j)} 

1]J — -L 



(1) 



Kenyon & Hartmann 1995 1 and do not take differences of the 
evolving surface gravities of PMS stars into account (private 
communication by the referee M.G. Hoare, and M.S. Bessel). somewhat equivalent to the usual \ i 



The sum is restricted to non-empty cells, and the nor- 
malisation by the simulated+observed density of stars in 
each cell gives a higher weight to more populated cellt^J 
Finally, the squared sum is divided by the total number of 

5 In Poisson statistics, where the uncertainty of a measure- 
ment M is a = VjV", our definition of Rrms turns out being 
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stars in the observed CMD, which makes R rms dimension- 
less, and preserves the number statistics when comparing 
clusters with unequal numbers of stars. 

The adopted strategy is: we start by building the Hess 
diagram corresponding to the J x (J — K s ) CMD of an 
actual young cluster. Next, we minimise Eq. [T] to find the 
values of P m in that give rise to the simulated Hess dia- 
gram that best matches the observation. Before reaching 
the absolute minimum, several ASA iterations are under- 
taken, which generates new starting points. At each new 
point (i.e, a new set of values P) produced by ASA, steps 
(i) to (xi) of Sect. [3] have to be repeated. 

By construction, the output of a single ASA run is the 
set of optimum values P m in (expected to correspond to 
the global minimum) and the respective R rm si but with 
no reference to uncertainties. However, given the statisti- 
cal nature of the process used to build the simulated Hess 
diagrams - as well as the stochasticity intrinsically associ- 
ated to young clusters, uncertainties to the parameters are 
strongly required. Thus, we run ASA several times (N run ) 
and compute the weighted mean of the parameters over all 
runs, using the R rms of each run as weight (w = 1/R rms ). 
It is important to remark that each new run begins with a 
totally new point P, with values randomly selected within 
the individual search-ranges (see above); this strategy pro- 
vides a means to minimise biases related to fixed initial con- 
ditions. Besides uncertainties, this procedure also provides 
a means to check for convergence patterns (see below). In 
summary, the ASA m i n approach consists of the full set of 
procedures so far described, namely, the star cluster sim- 
ulation (Sect. [3]), the N run ASA runs plus the statistical 
analysis (see above). Ideally, the parameters produced by 
ASA min should present a lower dispersion around the mean 
as Rrms declines. 
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Fig. 2. Convergence pattern for the models in Table[lJ with 
the input values shown by the vertical lines. Each symbol 
represents a single run. 



4.1. Tests with model star clusters 

We now employ template CMDs built with typical pa- 
rameters of PMS-rich young clusters to test the ability 
of ASA m i„ in recovering the input values. As described 
in Table [T] the models cover a broad range of proper- 
ties, with 150 < M du (M ) < 600, 5 < t c i u (Myr) < 25, 
3 < T sfs (Myr) < 16, < f bm < 1, including high 
DR absorption values (Model#5). In all cases we consider 
N S im — 100 and N run = 25 ASA runs. The retrieved pa- 
rameters, together with the corresponding R rmsi are listed 
in Table [TJ For simplicity, Table [l] lists only the param- 
eters found for the runs corresponding to the minimum, 
maximum and mean values of R rms . Note that the mean 
values are obtained by averaging out the N run independent 
outputs and using the respective R rms as weight; usually, 
the standard deviations are lower than 5% with respect to 
the mean values. The complete set of solutions is shown in 
Fig. [2} 

As expected, the retrieved parameters close in on the 
input values as R rm s declines, although with somewhat 
different dispersions around the mean among the param- 
eters. The convergence pattern is clearly tighter for M C \ U1 
E(J — K s ), and (m — M) j] followed by a(SAy) and t c i u , 
with fbin, T sfs, and 5 Ay at a third level. This is consis- 
tent with the mean values and respective standard devia- 
tions listed in Table [l] Nevertheless, it is remarkable that 
ASA m i n can retrieve the input values even for a CMD with 



such a high DR amount as that in Model#5. Also, its abil- 
ity to disentangle DR and binary fraction comes from the 
fact that, although both effects tend to redden the stellar 
sequences, binaries brighten them, while DR shifts them to- 
wards the opposite direction. It is also interesting the strong 
sensitivity for very short (e.g. Models#l and 4) and long 
(Models#2 and 5) SFSs. 

Similarly to Paper I, we also compute the cluster mass 
that potentially would be measured based on CMD prop- 
erties (Mcmd), which is an important piece of information 
for, e.g., establishing the dynamical state of a cluster. The 
first step is to find the actual mass of each star occurring 
in the CMD. However, this may be very difficult to find, 
especially in the presence of DR, SFS, unresolved binaries, 
and photometric uncertainties (which decreases the num- 
ber of stars that remain detectable in a CMD as the dis- 
tance modulus increases). In practise, what is usually done 
is: having derived the values of age, distance modulus and 
foreground reddening, Mcmd is computed by finding the 
probable mass for each star in the CMD by interpolation 
of the observed colour and magnitude among those of the 
nearest isochrones. Obviously, the precision of this proce- 
dure relies heavily on the amount of DR, photometric noise, 
binaries, etc. For instance, a heavily differentially reddened 
cluster of M c \ u ~ 600 M Q , characterised by a moderate bi- 
nary fraction and distance modulus, would have only ss 25% 
of its actual mass estimated based on CMD properties (e.g., 
Model#5 in Table [j}. 
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Table 1. Recovery of model parameters with ASA. 



Rank Rrms 


M c l u 


Age 


Tsfs 


(m - M)j 


E(J-K S ) 


SA V 


a(SA v ) 




McMD 




(M ) 


(Myr) 


(Myr) 


(mag) 


(mag) 


(mag) 


(mag) 




(M Q ) 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Model#l 


150 


25 


3 


9.0 


0.0 


0.5 


0.5 


0.0 




Min 0.1556 


149 


26.8 


3.7 


8.99 


0.00 


0.49 


0.51 


0.01 


133 


Max 0.1763 


155 


25.8 


3.8 


9.03 


0.01 


0.41 


0.55 


0.04 


145 


Mean 0.1653 


151 ±2 


25.9 ± 1.0 


3.3 ±0.6 


9.00 ±0.01 


0.01 ±0.01 


0.46 ± 0.05 


0.51 ±0.03 


0.01 ±0.01 


135 ±3 


Model#2 


250 


15 


13 


9.5 


0.5 


1.0 


1.5 


0.3 




Min 0.1124 


250 


15.3 


13.6 


9.50 


0.50 


0.94 


1.53 


0.30 


153 


Max 0.1240 


235 


15.6 


13.4 


9.45 


0.51 


0.70 


1.62 


0.25 


142 


Mean 0.1160 


250 ±6 


15.0 ±0.6 


12.9 ±0.7 


9.50 ±0.03 


0.50 ±0.01 


0.97 ±0.13 


1.53 ±0.05 


0.30 ±0.02 


153 ±5 


Model#3 


350 


10 


7 


10.0 


1.0 


1.7 


2.5 


0.6 




Min 0.1058 


354 


10.0 


6.5 


9.99 


1.00 


1.73 


2.54 


0.63 


200 


Max 0.1240 


355 


9.6 


4.8 


9.98 


1.00 


1.57 


2.87 


0.61 


204 


Mean 0.1105 


349 ±5 


9.9 ±0.5 


6.6 ±0.7 


10.00 ±0.03 


1.00 ±0.01 


1.61 ±0.16 


2.57 ±0.12 


0.60 ± 0.03 


193 ±3 


Model#4 


500 


5 


3 


10.5 


1.5 


2.5 


3.5 


1.0 




Min 0.1096 


503 


5.0 


3.0 


10.53 


1.51 


2.39 


3.59 


1.00 


278 


Max 0.1209 


503 


4.8 


2.9 


10.50 


1.49 


2.81 


3.25 


0.96 


261 


Mean 0.1143 


500 ± 11 


4.9 ±0.1 


2.9 ±0.3 


10.52 ±0.03 


1.50 ±0.01 


2.40 ±0.17 


3.53 ±0.11 


0.99 ±0.01 


280 ±6 


Model#5 


600 


20 


16 


10.3 


2.0 


10.0 


10.0 


0.8 




Min 0.1119 


609 


19.2 


15.9 


10.35 


2.01 


9.77 


10.28 


0.76 


154 


Max 0.1178 


569 


19.8 


16.1 


10.30 


2.00 


9.97 


8.80 


0.73 


140 


Mean 0.1149 


599 ± 14 


20.3 ±0.8 


16.3 ± 1.0 


10.31 ±0.02 


2.00 ±0.01 


9.93 ±0.26 


9.95 ± 0.40 


0.80 ±0.07 


152 ±4 



Table Notes. Cols. (1) and (2): R rms rank and value; Col. (3): actual cluster mass; Col. (4): cluster age; Col. (5): SFS timescale; 
Col. (6): apparent distance modulus in the J band; Col. (7): foreground reddening; Cols. (8) and (9): mean DR and standard 
deviation; Col. (10): unresolved binary fraction; Col. (11): cluster mass as measured in the CMD. The average stellar mass of 
the models is fnZ w 0.6 M . 



We close this section by concluding that ASA m i n - the 
minimisation of residuals between the observed and simu- 
lated Hess diagrams by means of ASA - is efficient in re- 
trieving the input parameters of model CMDs that cover a 
variety of conditions. 



5. Probing actual star clusters 

Based on the experience gained in Sect. |4.1| with model 
CMDs, we now move on to investigate actual yo ung clus- 
ters with ASA min . For this we use Collinder 197 jB onatto 



After several tests with ASA m i n , we settled on N S i m — 
50 and N run = 50 for Collinder 197 and N S i m — 500 and 



fc Bica||2010a[ ) and Pismis 5 ( |Bonatto fe Bica||2009 



BotTT 

clusters have been studied in Paper 1, which thus will al- 
low us to compare the effect of different assumptions on 
SFR, DR and binaries, on the derived parameters. They 
are projected near the Galactic equator (b ~ 1°) and have 
CMDs dominated by faint stars. To minimise confusion be- 
tween intrinsic PMS stars and red dwarfs of the Galactic 
field, we build field-star deco ntaminated CMD s by m eans 
of the algor ithm developed in Bonatto fc Bica (2007) and 
improved i n|Bonatto fc Bica ( 2010ap . As a result, the decon- 
taminated CMD of Collinder 197 has 690 stars (essentially 
PMS), while Pismis 5 has only 101, which is important for 
examining the effect of CMDs with different numbers of 
stars on our analysis. In addition, the colour- colour dia- 
gram s of both objects (Collind er 197: Fig. 7 of | Bonatto fc| 
Blca||2010a[ Pismis 5: Fig. 9 of T jBonatto fc Bica||2009bp "do 
not contain detections with abnormally high infrared ex- 
cesses that might be due to circumstellar material. 



25 for Pismis 5 as a compromise between conver- 



gence pattern and running timeQ Besides the normal DR 
distribution, for comparison purposes we also consider the 
case of a uniform (or flat) distribution. The parameters ob- 
tained with ASA m i n are given in Table [2] and the full set of 
runs is shown in Figs. [3] and [4] For comparison with other 
clusters, we also compute the bolometric magnitude and 
mass to light ratio. We remark that it is significant that 
convergence patterns reached after several thousand ASA 
iterations, similar to those of the model CMDs, occur for 
all parameters of both clusters. 

In both cases, the R rms corresponding to the normal DR 
distribution tends to be lower than those of the uniform dis- 
tribution. Also, reflecting the larger number of CMD stars, 
the R rms values of Collinder 197 are significantly lower than 
those of Pismis 5. 

Before moving on to interpreting the results, it is impor- 
tant to remind that each R rms value and corresponding op- 
timum parameters result from several thousand iterations, 
as ASA searches the R rm s hyper-surface for the absolute 
minimum of Eq. [T] At each iteration, N S i m twin clusters 
(i.e., consisting of exactly the same set of parameters) are 
built and incorporated into the simulated Hess diagram. 
Thus, the occurrence of a tight convergence pattern for the 



As a technical note we remark that each ASA run for 
Collinder 197 took ~ 1.2 hours, and ~ 2.5 hours for Pismis 5, 
on a single core of an Intel Core i7 920Q2.67 GHz processor. 
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Table 2. Parameters of Collindcr 197 and Pismis5 



lu A ge r sfs (m-M)j E(J - K s ) SA V a(6Ay) f bin M C md d Q 10 3 MLR M BOL 

(M Q ) (Myr) (Myr) (mag) (mag) (mag) (mag) (M ) (kpc) (M /L Q ) (mag) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

Collinder 197 - Differential reddening mode: normal 

0.4376 381 12.9 12.1 8.78 0.01 2.55 1.24 0.00 211 0.57 4.7 -7.5 

0.4442 344 12.0 11.5 8.82 0.00 2.63 1.08 0.02 182 0.58 6.6 -7.0 

0.4401 376 ±13 13.6 ± 0.9 13.1 ± 1.0 8.78 ±0.07 0.03 ± 0.02 2.44 ± 0.15 1.20 ± 0.06 0.01 ± 0.01 208 ± 7 0.55 ± 0.02 4.2 ± 1.0 -7.6 ±0.2 

Alrr ^7fi+ 230 12 9+ 7 ' 2 121+ ' 7 8 78+ 1 - 30 001+ 017 2 55+ 1 ' 02 1 24+ ' 96 00+ ' 39 211 + 30 57+ - 33 A 7+ 3 - 3 _7 ^0.1 

^'D-igO lz - y _4.7 lz - i -5.3 s -' 5 -0.76 u - ui -0.01 z ' oa -0.96 J " z *-0.51 u - uu -0.00 Zii -56 u - a '-0.17 4 -'-0.4 1 -°+0.6 

la 396 ±96 14.3 ±3.5 10.0 ± 1.7 8.95 ± 0.47 0.01 ± 0.01 2.58 ± 0.50 1.42 ± 0.40 0.17 ± 0.11 224 ± 12 0.58 ± 0.12 3.5 ± 1.9 -7.9 ±0.6 

Collinder 197 - Differential reddening mode: uniform 

0.4420 375 13.8 13.6 8.94 0.12 3.87 0.01 225 0.56 5.3 -7.4 

0.4501 364 15.5 15.0 8.72 0.10 3.84 0.00 224 0.51 4.7 -7.5 

0.4453 369 ± 12 13.8 ± 1.1 13.5 ± 1.1 8.94 ±0.09 0.13 ± 0.01 3.76 ± 0.11 — 0.01 ± 0.01 227 ± 8 0.55 ± 0.02 4.1 ±0.8 -7.6 ±0.2 

Pismis 5 - Differential reddening mode: normal 

0.5674 145 7.3 5.3 10.48 0.14 0.63 2.25 0.75 80 1.11 3.1 -6.9 

0.5729 147 4.0 1.2 10.77 0.11 0.44 2.80 0.83 100 1.30 1.8 -7.5 

0.5692 142 ±7 6.8 ± 0.8 5.0 ±1.0 10.56 ± 0.08 0.13 ± 0.01 0.56 ± 0.29 2.30 ± 0.25 0.83 ± 0.08 81 ±4 1.16 ± 0.04 2.6 ± 0.4 -7.1 ±0.2 

Ala 14 5 ±H 7.31*;° 5.3±^ 7 10.481-°;?! 0.14+°;°? 0.63t£;» 2.25^ 0.75+°$ 80+J 23 l.lllg^ 3.11*3 

la 154 ± 50 9.0 ±3.2 5.0 ±1.3 10.52 ± 0.38 0.14 ± 0.04 0.92 ± 0.54 2.62 ± 0.70 0.73 ± 0.16 82 ± 7 1.12 ± 0.20 2.7 ±1.5 -7.0 ± 0.6 



Pismis 5 - Differential reddening mode: uniform 



0.5696 110 12.2 11.5 9.82 0.13 3.85 0.01 57 0.83 4.5 -6.2 

0.5743 125 12.2 11.2 10.01 0.16 3.89 0.02 63 0.88 4.5 -6.4 

0.5726 110 ±8 13.0 ± 1.1 12.2 ± 1.2 9.84 ±0.08 0.14 ± 0.01 3.79 ± 0.13 — 0.05 ± 0.04 57 ± 4 0.83 ± 0.03 4.4 ± 0.9 -6.3 ± 0.2 
Table Notes. Col. (1): minimum, maximum and mean R rm s values; Col. (2): actual cluster mass; Col. (3): cluster age; Col. (4): SFS; Col. (5): 
apparent distance modulus in the J band; Col. (6): foreground reddening; Cols. (7) and (8): differential reddening; Col. (9): unresolved 
binary fraction; Col. (10): cluster mass as measured in the CMD; Col. (11): distance from the Sun; Col. (12): bolometric mass to light 
ratio; Col. (13): bolometric magnitude. Alcr: parameters occurring at the boundaries of the lcr domain, lcr: weighted average of values 
occurring within the la domain. 



absolute minimum parameters over a series of independent 
runs shows a self consistency of the method and cannot be 
taken as fortuitous or model dependent. Instead, it would 
be strongly indicative that the optimum parameters pro- 
duced by ASA m i n are indeed representative of those of 
the cluster being studied. On the other hand, we remark 
that this argument applies only to the absolute minimum 
of each ASA m i n run, since the R rm s topology around this 
feature is not taken into account. In this sense, the quoted 
errors should be taken as internal, probably not reflecting 
the realistic parameter uncertainties. Indeed, as we show 
in Sect. |5.2| in some cases the drop towards the absolute 
minimum is somewhat gentle, which means that there can 
be a significant dispersion (different values at similar R rms 
levels) around the absolute minima. Although the depres- 
sion shape departs from gaussianity, we estimate the ap- 
proximate extension of the lcr domain and compute the 
weighted mean and dispersion for the parameters occur- 
ring inside it (again using the individual R rm s as weight). 
The results - restricted to the uniform DR - are given in 
the additional entry labelled as lcr in Table [2| Compared to 
the previous statistics, the uncertainties now are more re- 
alistic and consistent both with the R rm s topology and the 
assumptions incorporated into the simulations. For com- 
pleteness, we also provide in Table [2] the parameter values 
that occur at the boundaries of the lcr domain (labelled 
as Alcr). In general, they imply considerably larger errors 



than before. However, given the underlying assumption of 
residual minimisation and respective R rm s value of different 
parameter sets, we believe that the weighted-mean values 
(la) are more representative for the clusters dealt with in 
this work. In addition, the mean values arc essentially un- 
changed, except for a somewhat higher binary fraction for 
Collinder 197. 

Collinder 197: The values of cluster mass, age, SFS, bi- 
nary fraction, and distance from the Sun are essentially 
insensitive to the DR mode. They consistently indicate a 
cluster of M clu « 400 M Q , t clu a 14 Myr, r sfs « 10 Myr, 
with a relatively low binary fraction of fu n is 0.2, located 
at ~ 0.6 kpc from the Sun. However, assuming that repre- 
sentativeness increases for lower R rms values, DR in this 
cluster would follow a normal distribution characterised by 
SAy « 2.6 and a (5 Ay) ~ 1.4. The differences that occur 
in distance modulus and foreground reddening can be ac- 
counted for by the single, fixed value of 5 Ay of the uniform 
DR distribution, which requires higher values of (m — M)j 
and E(J — K s ) than the normal mode to account for the 
observed colour spread of the stellar sequences. The fore- 
ground reddening is very low, corresponding to Ay w 0.1. 

Pismis 5 : Compared to Collinder 197, here the DR modes 
result in significant differences for most of the cluster pa- 
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Fig. 3. Similar to Fig. |gj for Collinder 197. Results for 
the uniform (empty symbols) and normal (filled) differ- 
ential reddening modes are shown. Vertical lines corre- 
spond to the mean values for the uniform (dotted) and 
normal (solid) distributions. The abscissa scales represent 
the search-ranges used at the beginning of the ASA min- 
imisation process. 



rameters, except for the foreground reddening. Also, as ex- 
pected from the lower number of stars, the convergence pat- 
terns are in general looser than for Collinder 197. Perhaps 
the most contrasting result lies in the binary fraction that 
reaches fun ~ 0.7. Lacking one degree of freedom, the uni- 
form DR mode tries to describe the CMD spread by as- 
suming a large value of 8 Ay with essentially no binaries. 
In contrast, the normal mode requires a lower 8 Ay with 
a high dispersion a (8 Ay), which yields a high binary frac- 
tion. Combined, the low 8 Ay and high /;,,„ of the normal 
mode put Pismis 5 at a distance ~ 40% larger than that 
implied by the uniform mode. The foreground reddening, 
Ay 0.9, is higher than that of Collinder 197. 

Interestingly, because of the additional free parameter, 
the normal DR distribution produces a cluster age signifi- 
cantly younger than that implied by the uniform distribu- 
tion, especially for Pismis 5. Also, in both cases the SFS is 
equivalent to ~ 60% of the cluster age. 

Comparing the values obtained by ASA min with those 
in Paper I, we see that now the ages tend to be some- 
what older and the distances shorter, with a reasonable 
agreement among the other parameters. A similar conclu- 
sion applies to the scarce works on both objects. Previous 
estimates for Collinder 197 are: age~ 5 ± 4 Myr, d® ~ 
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Fig. 4. Same as Fig. [3] for Pismis 5. 



1.1 ± 0.2 kpc, and mass 100 - 500 M jBonatto & Bica 
2010a and references therein). For Pismis 5, they are: age 
5 - 15 Myr and ri - 1.0 ±0.2 kpc l ]Bonatto fc Bica||2009b 
and references therein) . In paper I we had assumed a con- 
stant SFR, uniform DR and a fixed binary fraction. Despite 
the latter two effects, the main source of differences is the 
SFR. In a constant SFR, stars of any age have the same 
probability of being formed, and young stars are signifi- 
cantly brighter than their older counterparts. Thus, when 
trying to match the same observed CMD, a simulation that 
contains an enhanced population of young stars (constant 
SFR) would require higher values of distance modulus than 
another based on a decaying SFR. At the same time, the 
constant SFR simulation would also imply a younger age. 

The representativeness reached by the minimisation 
process described above for Collinder 197 and Pismis 5 can 
be visually appreciated by comparing the observed and 
simulated Hess diagrams (Fig. [5]). The latter have been 
constructed with the mean parameters and assuming the 
normally-distributed DR (Table [2). Although the occur- 
rence of some discreteness in both Hess diagrams, which is 
a natural consequence of the relatively low-mass nature of 
the clusters, especially for Pismis 5, observation and simu- 
lation show a reasonable correspondence in both clusters. 

Finally, in Fig. [6] we compare a single CMD realisation 
randomly selected among the N sim simulated clusters - but 
having exactly the same number of stars - with the observed 
CMDs. The DR distributions (Collinder 197 has higher My 
and lower a(SAy) than Pismis 5) are also shown. For il- 
lustrative purposes, Fig. [6] also shows the isochrones that 
represent the full star-formation history of both clusters. 
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Fig. 5. Observed (top panels) and simulated (bottom) 
Hess diagrams of Collinder 197 (left) and Pismis5 (right). 
Lighter grey shades indicate higher densities of stars. The 
simulated diagrams have been built with the mean param- 
eters and the normal DR mode (Table [2]) . 



Fig. 6. Field-star decontaminated CMDs (light-shaded cir- 
cles) of Collinder 197 (left) and Pismis 5 (right) com- 
pared to a single realisation taken from the respective 
simulations (empty circles). The extraction radii are 10' 
(Collinder 197) and 6' (Pismis 5). Isochrones representing 
the star-formation history have been set according to the 
mean parameters (Table pi). Insets show the DR distribu- 
tion. 



They have been set with the mean values of distance mod- 
ulus and foreground reddening in Table [2j Note that the 
optimum simulation and respective isochrone solution end 
up naturally respecting the blue border of the observed 
stellar sequences as a (not-imposed) boundary conditiorQ 
Also, the fading and reddening effect of DR on the stars 
is clearly seen when one compares the youngest (and red- 
dest) isochrone with the redwards spread of the PMS stel- 
lar sequences. This also shows that, if DR is not properly 
taken into account, fitting isochrones to a CMD would re- 
quire somewhat younger ages coupled to higher values of 
distance modulus and foreground reddening, especially to 
account for the faint and red PMS stars together with the 
blue border. 

Again, given that both are low-mass, young clusters, 
some morphological differences should be statistically ex- 
pected, especially in the MS. Nevertheless, simulated and 
observed CMDs are similar in both cases. 



7 From our perspective, this is somehow reassuring, since the 
blue border has been taken as a constraint to estimate funda- 
mental pa rameters of young cluste rs with simpler methods such 
as that in Bonatto & Bica (2010a I. 



5.1. The low-mass stochasticity 

An interesting issue that can be investigated with ASA m i n 
is the natural stochasticity associated with low-mass clus- 
ters. In other words, how the number of simulated clusters 
{N 8 i m ) affect the convergence pattern of the retrieved pa- 
rameters. Obviously, N sim is expected to play an important 
role especially on a poorly-populated (M c i u < 150 M ) and 
young (t c i u r~j 7Myr) cluster such as Pismis 5, in which the 
stochasticity tends to be critical. 

For a low-mass cluster consisting essentially of PMS 
stars, the stochasticity issue can be summarised as follows. 
Statistically, a random simulation of this cluster (same age, 
mass, etc) may contain a massive, bright star that is not 
present in the actual cluster. Consequently, because of the 
mass constraint, the CMD of this particular simulation 
would also lack a large fraction of the low-mass content. 
Thus, despite having been built with exactly the same pa- 
rameters as the cluster, the R rms residuals of this simula- 
tion would be large, with a low representativeness. Given 
the extremely large number of parameter combinations, the 
probability of finding a single simulation with a Hess dia- 
gram matching that of the cluster is vanishingly low. On 
the other hand, when many simulations with exactly the 
same parameters are considered, the individual Hess dia- 
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Fig. 7. Increasing the number of twin simulations (N s i m ) 
improves the convergence pattern for a poorly-populated 
cluster such as Pismis5. Simulations built assuming the 
normal DR mode. The number of simulations used in the 
runs are N S i m = 10 (light-shaded circles), N S i m = 100 
(empty), and N sim — 500 (filled). 

grams would certainly present different representativeness. 
Thus, the average Hess diagram over many simulations is 
expected to present a higher similarity (i.e., low R rm s) with 
that of a cluster's than what would be obtained with a sin- 
gle, random simulation. It is in this context that averaging 
out several simulations of such twin clusters becomes an 
effective way of minimising stochasticity. 

Such effect is illustrated in Fig. [7J in which we show 
the retrieved parameters for Pismis 5 after running ASA m i n 
with N S i m = 10, 100 and 500, always assuming the same 
initial conditions. Clearly, the convergence pattern, which 
is very loose for N sim = 10, becomes tighter as N sim in- 
creases. Despite the significant scatter associated mainly 
with N S i m = 10, it is interesting to note that at the lowest 
Rrms values, the parameters of N S i m = 10 and 100 converge 
to those of N S i m = 500. Besides, the differences between the 
N s im — 100 and N S i m = 500 runs are significantly smaller 
than those between N S i rn = 10 and N s i m = 100. 

Note, however, that some scatter is still present for SFS, 
mean DR and dispersion, and binary fraction, even when 
Nsim = 500. Ideally, a very large N sim should be used to 
reach a very tight convergence pattern for a low-mass clus- 
ter. In practice, however, this may lead to exceedingly long 
runtimes. In this sense, our recipe of using a moderate N sim 
combined to a series of runs may be taken as a compromise 
between runtime and robustness. 



Fig. 8. Selected two-dimensional R rms projections for 
Collinder 197 (left) and Pismis 5 (right). Darker colours in- 
dicate lower Rrms- Convergence patterns vary significantly 
among different parameters. 

5.2. ASA m ,in and the R rm s topology 

Having found the best-fit parameters, we now use them to 
examine the topology of the R rms hyper-surface, something 
we address by means of selected two-dimensional projec- 
tions (Fig. pjb. For statistical significance, the maps were 
produced with N s i m = 100 and 500 and the absolute min- 
imum R r ms parameters, respectively for Collinder 197 and 
Pismis 5. The presence of a minimum is clear, but with con- 
vergence patterns varying significantly among the parame- 
ters, being tight for most but somewhat loose especially for 
the binary fraction. 

Interestingly, the projections now are significantly 
smoother and with less features than the equivalent ones 
shown in Fig. 5 of Paper I. Possible reasons for such a con- 
trast are that in Paper I the R rms hyper-surface (i) was 
built using some fixed parameters (binary fraction and SFS) 
and restrictive conditions (flat SFR and uniform DR distri- 
bution), (ii) the free parameters were allowed to vary within 
bins of fixed size (and not continuously as in the present 
approach) and within restricted ranges, and (Hi) we used a 
different statistics for finding solutions. Apparently, the ad- 
ditional free parameters and less-restrictive conditions have 
raised the degeneracy of solutions found in Paper I. 

In any case, what is really important for our approach 
is the existence of at least one conspicuous minimum, and 
our expectation (Sect. [4]) is that ASA is efficient in find- 
ing its way towards the absolute minimum. In addition, 
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Fig. 9. Selected one-dimensional projections of R rm s com- 
puted with N S i m — 500 for Collinder 197 (left) and N S i m — 
1500 for Pismis5 (right). Besides the values computed for 
the absolute minimum parameters (solid line) , we also show 
the Rrms projections using values 10% lower (dotted) and 
higher (dashed) than the minimum. The absolute minimum 
is shown by the shaded circle, and the approximate la do- 
main is shown by the grey region. 



a single minimum also serves to strengthen the unicity of 
the solutions. Thus, given the relevance of this assumption 
to the results discussed above, we summarise it in Fig. [9] 
with one-dimensional projections of the R rms hyper-surface 
for selected parameters of both clusters. For instance, we 
take from Table the minimum values (obtained with the 
normal DR mode) for all the parameters except cluster 
mass, which is allowed to vary within a wide range (in this 
case, 0.5 < M c [ u (Mq) < 250). Then, we compute R rms for 
masses within the adopted range, but keeping the remain- 
ing parameters hxed. For a deeper perspective on the effect 
of parameter variation on the R rm s shape, this step is re- 
peated with the fixed parameters changed to 10% higher 
and lower than the optimum values. The same procedure is 
applied to the age, mean DR and binary fraction. 

Consistently with the R rms maps (Fig. [8]), the selected 
Rrms projections present a single minimum with a degree of 
definition (depth and width) that varies significantly among 
the adopted thresholds. When comparing different param- 
eters, the minima tend to be quite narrow and deep for the 
age, somewhat wide and shallow for the binary fraction, 
and intermediate for the remaining parameters. 



The general conclusions emerging from the above dis- 
cussion are: (i) the R rm s hyper-surface contains at least one 
minimum within the adopted search range, (it) the morpho- 
logical features of the minimum vary according to each pa- 
rameter, and most importantly, (Hi) ASA m i n finds its way 
through the i?, rms topology towards the absolute minimum. 



6. Summary and conclusions 

A new approach (ASA m i n ), designed to obtain a set of im- 
portant parameters of young star clusters at a statistically 
significant confidence level, is presented in this paper. In 
short, it is essentially based on photometric properties and 
involves building realistic simulations of the Hess diagram 
of an actual star cluster, from which the residuals {Rrms) 
with respect to the observed Hess diagram are computed. 
Besides cluster mass, age, foreground reddening and dis- 
tance modulus, the simulations include the SFS, (alterna- 
tive modes of) DR and the binary fraction as free param- 
eters. The CMD spread due to photometric uncertainties 
is explicitly taken into account. Important features of the 
simulations are a linearly decreasing SFR and a normally 
distributed DR. 

To find the absolute minimum of the R rm s hyper-surface 
we use the global optimisation method known as adap- 
tive simulated annealing (ASA), which is rather efficient 
and capable of escaping from local depressions. Given the 
highly-statistical nature of the simulations (and the clus- 
ter stochasticity) , we show that an acceptable parameter- 
retrieval rate is achieved by combining a moderate number 
of simulated clusters with a series of independent ASA runs, 
while realistic errors in the derived parameters are obtained 
by exploring properties of the depression shape. Tests with 
model clusters built with a broad range of parameters show 
that the distribution of retrieved values (corresponding to 
the absolute R rm s minimum) usually follows a convergence 
pattern that is tighter as R rm s declines, with the same oc- 
curring with actual clusters. We also find that the param- 
eter retrieval presents a high sensitivity for cluster mass, 
distance modulus and foreground reddening, but dropping 
somewhat for the remaining parameters. 

We remark that the particular results discussed in this 
work may be somewhat model dependent, in the sense that 
they are based on the 2MASS near-infrared photometry 
coupled to the Padova and Siess MS and PMS isochrone 
sets. Other isochrones with somewhat different mass to light 
ratios for individual stars may possibly affect the star clus- 
ter parameters when retrieved by ASA m i n , especially the 
age, mass and distance. In any case, ASA m in can be eas- 
ily adapted to any photometric system (and isochrone set), 
provided the respective MS and PMS isochrones and the re- 
lation of photometric errors with magnitude are available. 
However, by allowing a deeper view through the dust, the 
near-infrared seems to be the best window to disentangle 
the effects of DR, binaries and SFS. 

In summary, we show in this work that our semi- 
analytical and comprehensive approach is capable of uncov- 
ering a series of parameters of young clusters, even when 
photometry is the only available information. Among these, 
the cluster age, star-formation spread, mass and binary 
fraction are important for establishing the dynamical state 
of a cluster, or derive a more precise SFR in the Galaxy. 
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